Partial order in Potts models on the generalized decorated square lattice 
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We explore the Potts model on the generalized decorated square lattice, with both nearest (Ji) 
and next-neighbor (J2) interactions. Using the tensor renormalization-group method augmented 
by higher-order singular value decompositions, we calculate the spontaneous magnetization of the 
Potts model with q — 2, 3, and 4. The results for q = 2 allow us to benchmark our numerics using 
the exact solution. For g = 3, we find a highly degenerate ground state with partial order on a single 
sublattice, but with vanishing entropy per site, and we obtain the phase diagram as a function of 
the ratio J2/J1. There is no finite-temperature transition for the q = 4 case when Ji — J2, but the 
magnetic susceptibility diverges as the temperature goes to zero, showing that the model is critical 
at T = 0. 

PACS numbers: 64.60. Cn, 05.50.-|-q, 75.10.Hk, 64.60.F- 



The (j-state Potts model plays an important role in the 
understanding of different phases and critical phenomena 
(for a review see Ref. [|l[). The behavior of ferromagnetic 
Potts models is well understood due to their universal- 
ity. The properties of antiferromagnetic Potts models are 
much more complex, varying widely for different values 
of q and lattice topologies. On every lattice there exists a 
value Qc such that for q > Qc the model has exponentially 
decaying correlations and hence no order at all tempera- 
tures, including T = 0. For q = qc, the model has a crit- 
ical point at zero temperature, and for q < q^ all types 
of behavior are possible Q , including long-ranged order 
or partial order. When the lattice is irregular, meaning 
that not all sites are equivalent, entropy-driven phase 
transtions and exotic forms of partial order may arise 
P-H . Partial order has been found both in classical sta- 
tistical models p, M] and in quantum systems H . 

The decorated square lattice is formed by adding one 
site to each bond of a square lattice. Here we study 
Potts models on the generalized decorated square lattice 
obtained by introducing an interaction between the dec- 
orating sites [Fig. 0(a)]. This is a tripartite lattice, for 
which one attribution of sublattices {A, B, C) is shown 
in Fig. H(b). The Hamiltonian of the model is 





H 






({i.j)) i(^C 



where ai = 0,1,2,... is the Potts state on lattice site 
i. The first term denotes all nearest- neighbor bonds and 
the second all next- neighbor bonds in Fig. 0(a), while the 
third term denotes a field term h applied to the sites in 
any sublattice(s) C for the calculation of magnetizations 
and susceptibilities. 

In the pure decorated square lattice, with only nearest- 
neighbor ( Ji) interactions, taking a partial trace over the 
decorating {B and C) sites effects a mapping to a ferro- 
magnetic Potts model on the square lattice ]9|-pd| . From 
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FIG. 1: (a) Generalized decorated square lattice with interac- 
tions between decorating sites, (b) Sublattice labels for this 
tripartite lattice, (c) Label definition for the tensor expressed 
in Eq. (|). 



this one may deduce that the critical value for the g-state 
antiferromagnetic Potts model on the decorated square 
lattice is qc = \{?i + ^/b) = 2.618, meaning that for q = 3 
and q — A there is no ordered phase at low tempera- 
tures and no phase transition occurs. This value of qc is 
smaller than that for the antiferromagnetic Potts model 
on the square lattice, where qc = i p^-p^. For q = 2 
(the Ising model), the exact solution is one with complete 
long-range order [ p7| . 

When next-neighbor interactions are introduced, reen- 
trant behavior of the ordered phases has been found both 
for g = 2 [|l|, |l|l and g = 3 (l|l. However, this phe- 
nomenon required that one of the interactions be ferro- 
magnetic Ip^ . Here we consider the regime where both 
Ji and J2 are antiferromagnetic and focus primarily on 
the isotropic situation (Ji = J2), where the competition 
between the different interactions is the largest. 

For a classical model on any lattice, it is always pos- 
sible to find a tensor-network representation of the par- 
tition function pG, Gl|. For the Hamiltonian defined in 
Eq. (m), we define a rank- four tensor Tijki for the unit 



cell shown in Fig. |l|(c) as 

Tijkl = ^CXp[-/3Jl((5cr,,o 

m 

exp[-/3J2((5o-.,o 



„+ <^<Tfc,CT„+(5cr,,cr„)]x 

. + Sa,,^, + 6„,,^J]. (2) 



The Potts variables (Ji on each site serve as the indices 
of the tensors, and the tensors for each unit cell form a 
square-lattice tensor network. 

Recently, we have developed a Tensor Renormaliza- 
tion Group (TRG) technique which makes use of higher- 
order singular value decomposition (HOSVD) of the ten- 
sors to truncate them during the renormalization (coarse- 
graining) process |2^ . There are two steps in this method 
at each iteration n. In the first, we block two tensors to 
form a larger tensor. 
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where a; = xi (^ a;2 and x' = x'l (g) Xg. In the second, we 
perform a HOSVD, 

M^^Xv = E S^JklU^.U^XkU!yl, (4) 

ijkl 

where the U matrices are unitary and S is the core tensor 
of M("\ which is fully orthogonal and pseudo-diagonal 
|P3| . The tensor is then truncated according to 



T, 



(n+l) 
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to obtain a new tensor r("+i) with the same dimension 
as r("). In one "HOTRG" step, the size of the lattice 
is reduced by a factor of 1/2. By repeating the HOTRG 
step alternately in the x and y directions, we contract 
the tensor network and obtain the partition function of 
the model. 

While several methods exist to contract the tensor net- 
work [go, 21, 2^, ^, we have demonstrated |2^ that the 
accuracy of this HOTRG method is very high: a tensor 
dimension of just D = 24 in the two-dimensional Ising 
model is sufficient to ensure that the relative error in the 
free energy is less than 10""^, even at the critical tempera- 
ture. The accuracy can be further improved by including 
the effects of the environment in the renormalization pro- 
cess, known as the second renormalization group (SRG) 
method p4| , although we do not need to implement this 
for the accuracies attained here. 

For numerical reasons, the thermodynamic quantities 
obtained most accurately from the HOTRG partition 
function are the magnetization and susceptibility. The 
small field term included in Eq. (|f|) pins a state oi q = 
order to a specific sublattice £. For a starting square 
lattice of dimensions LxL = N, we define the order pa- 
rameter as 



^^ = ^(E^-o)-^' 



(6) 
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FIG. 2: Magnetization of the q = 2 Potts model on the gen- 
eralized decorated square lattice with Ji = J2 = 1. The bond 
dimension in the HOTRG calculation is D = 30. An Ising-like 
phase transition occurs at Tc = 0.97423 



which is the expectation value of the number of spins 
in sublattice L whose value is 0. The term \/q is sub- 
tracted to ensure a value of in the disordered phase. 
The susceptibility is given by x = dM/dh. The order 
parameter M expresses the broken sublattice symmetry 
p2[ of the partially ordered phase: only sites of one sub- 
lattice are ordered, while all others are disordered (but 
may still be correlated locally). An ordered state with 
broken sublattice symmetry was first suggested for the 
q = 3 antiferromagnetic Potts model on the square lat- 
tice, but more detailed analysis demonstrated that this 
order is artificial, the model being disordered at all tem- 
peratures other than a Tc = critical point |I2-n3]. 

The generalized decorated square lattice model is ex- 
actly solvable |I8 for q = 2. There are three sublat- 
tices but only two possible states for each site, leading 
to frustration. Because the coordination number of the 
A sites, ZA ~ 4, is less than that of the B and C sites 
{zb = zc = 6), the ground state is one with Neel order 
on the B and C sublattices while the Ising spins on the 
A sites are free. This configuration minimizes the num- 
ber of unsatisfied bonds and hence the energy. The free 
spins are responsible for the zero-temperature residual 
entropy, Sq = | In 2 per site. The universality class of 
the phase transition from the Neel phase to the param- 
agnetic phase remains that of the two-dimensional Ising 
model. In Fig. we show the magnetization obtained for 
Ji = J2 with tensor dimension D = 30 in the HOTRG 
calculation. The transition temperature we compute is 
Tc = 0.97423, which is very close to the exact result 
Tc ~ 0.9742197, demonstrating clearly the accuracy of 
the HOTRG method. 
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FIG. 3: Ground-state configurations for the q — 3 Potts model 
on the generahzed decorated square lattice. The three states 
are represented by red for 0, black for 1, and green for 2. Pan- 
els (a) and (b) illustrate an arbitrarily chosen square on the B 
and C sublattices with one pair of diagonally opposite sites in 
the same state, which creates a complete partial order on the 
corresponding sublattice. In panel (c), both opposite pairs 
have the same state and all three sublattice are ordered. All 
of the ground-state configurations are summarized in panel 
(d) and a n/2 rotation of this. 



We turn now to the g = 3 model and begin by develop- 
ing a heuristic understanding of the underlying physics. 
Because each site may take one of three states and the 
lattice is tripartite, all bonds can be satisfied and this 
model is unfrustrated. If the three states are represented 
by three different colors, as in Fig. |^, the problem of find- 
ing all the ground-state configurations may be cast as a 
coloring problem. In fact the division of sites into three 
sublattices shovifn in Fig. |l|(b) is not unique (a situation 
that would ensure the same type of order as in the q = 3 
Potts model on the triangular lattice). The origin of the 
different possible sublattice designations is the fact that 
Za = 4: for one of the sublattices, introducing more de- 
grees of freedom and causing a degenerate ground state. 

For the square of B- and C-sublattice sites marked 
by the orange dashed line in Figs. ||(a) and (b), at least 
one pair of diagonally opposite sites must have the same 
color. Each such square effectively nucleates a state of 
complete order on the relevant sublattice, as shown in 
Figs. ||(a) and (b), without ordering the other two sub- 
lattices (each of which in fact forms a state of two- fold de- 
generate alternating lines). If both pairs of opposite sites 
have the same color [Fig. ^(c)], then the nucleated state 
would be one of complete order on all three sublattices, 
but this is a set of measure zero. The set of all ground- 
state configurations is summarized in Fig. 0(d), and it is 



one with partial order on a single sublattice. Thus the 
additional freedom contained in the special topology of 
the generalized decorated square lattice is not sufficient 
to drive the system into a disordered phase. 

To investigate the nature of the partially ordered state, 
if the B sublattice is ordered with g = 0, the state 
of the intervening lines of A and C sites may be ei- 
ther 1212 ... 12 or 2121 ... 21 at random, and similarly 
for an ordered C sublattice. Thus the partially or- 
dered state breaks the 7r/2 rotation symmetry of the lat- 
tice. The ground-state degeneracy at zero temperature 
is iVo = 6(2^ — 1), where 6 is subtracted because config- 
urations in which both B and C sublattices are ordered 
[Fig. ^(c)] are double-counted. Despite this high degner- 
acy, the residual entropy per site. 



5'o = lini 



log(6 X 2^ - 6) 
L2 



= 0, 



(7) 



vanishes because of the one-dimensional nature of the 
freedom in ground-state configurations within a two- 
dimensional lattice. This type of partial order is dif- 
ferent from the order found in Refs. B, p2l, where only 
one sublattice was ordered and the other sites remained 
in random states. For the generalized decorated square 
lattice, the sites on lines between the ordered sublattice 
sites remain highly correlated, forming a one-dimensional 
antiferromagnetic Ising model {q = 2), which has two- 
fold degeneracy between the alternating states as noted 
above. 

We calculate the magnetization and susceptibility for 
the g = 3 Potts model on the generalized decorated 
square lattice for Ji = J2 using the HOTRG method 
with D = 30. As shown in Fig. ^, we find that the mag- 
netization goes continuously to zero at a second-order 
phase transition at Tc = 0.528(1). At low temperatures, 
the magnetization defined by Eq. (^) with C = C takes 
the value 2/3, which means that all C sites are fully po- 
larized and confirms our analysis concerning the nature 
of the ground state. The susceptibility (inset. Fig. Q) 
shows a clear divergence at Tc, confirming the second- 
order nature of the phase transition. 

Having established the physics of partial order where 
Ji = J2, we next investigate the regime where the bonds 
are not unequal (Ji ^0/2)- We consider only the antifer- 
romagnetic region where Ji, J2 > 0. In the limit where 
J2 = 0, the system is a pure decorated square lattice with 
qc = 2.618 and a disordered ground state for q ~ 3. In 
the limit where Ji = 0, the system is a g = 3 Potts model 
on the square lattice formed by the decorating B and C 
sites, while the A sites become isolated. For the antifer- 
romagnetic Potts model on the square lattice, the critical 
value is gc = 3 [^, and hence Tc = here. The phase 
diagram we calculate is presented in Fig. pi with the crit- 
ical temperature shown in units of Ji when Ji > J2 and 
J2 when J2 > Ji- 
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FIG. 4: Magnetization for the q = 3 Potts model on the 
generalized decorated square lattice with Ji — J2 and D = 30. 
A second-order phase transition occurs at Tc — 0.528(1). The 
inset shows the susceptibility, which diverges at Tc. 



FIG. 5: Phase diagram for the q — 3 Potts model on the 
generalized decorated square lattice as a function of J2/J1. 
HOTRG calculations for the critical temperature Tc were per- 
formed with D = 30. The temperature is expressed in units 
of Ji for J2/J1 < 1 and in units of J2 when J2/J1 > 1. 



Considering first the situation J2 < Ji, we find a finite- 
temperature phase transition to a state of partial order 
for all J2, and that Tc — > continuously as J2 — > 0. This 
implies that Qc > S in this limit. However, it is known 
that Qc = i(3 + V5) = 2.618 at J2 = [|-0, mean- 
ing a disordered ground state at T = 0, and hence we 
find discontinuous behavior here. The key property is a 
qualitative one, the connectivity of the lattice, and not a 
quantitative one related to the coupling ratio. This result 
may in fact be anticipated from the fact that the ground 
state should be the same in the whole region Ji, J2 > 0, 
because it corresponds to the same coloring problem, and 
hence Tc should be finite for all J2 > 0. 

We have not been able to find a similar example of 
such a discontinuity anywhere in the literature. We note 
in this context that the applied field acts to induce or- 
der at finite temperatures in the antiferromagnetic Ising 
model on the triangular lattice and the g = 3 Potts model 
on square lattice [g7| , both of which are critical at T = 
in the absence of a field. In principle our method is sus- 
ceptible to this type of physics, because our calculations 
use a small applied field, and so we have verified our re- 
sults (Fig. m by varying h. Furthermore, our system has 
the essential difference that the ground state of the deco- 
rated square lattice is not critical but is fully disordered 
(exponentially decaying correlations) in the absence of a 
field. At the other side of the phase diagram in Fig. ||, 
obtained by taking Ji —>■ 0, the model becomes critical 
at zero temperature and fiuctuations become very large 
WM. In Fig. H we show the critical temperature calcu- 
lated for moderate ratios J2/J1, and we appeal to known 
results in expecting that Tc continues smoothly to zero 



when Ji — > 0. 

We conclude our analysis of the generalized decorated 
square lattice with the q = A case, where we consider only 
J2 = Ji- For this value of 9, the antiferroiriagnetic Potts 
model on the Union- Jack lattice was found to have a 
finite-temperature phase transition to a partially ordered 
ground state |Q, g]. The generalized decorated square 
lattice with J2 = Ji may be considered as the result 
of eliminating half of the fourfold-coordinated centering 
sites within the squares of the Union- Jack lattice. This 
weakens the constraints on the ground state, which thus 
has sufficient freedom that a disordered state may well 
be expected: there are four states per site and the lattice 
is tripartite, with many ways to divide the system into 
three sublattices. 

We compute the magnetization by the HOTRG 
method, finding no order over the whole temperature 
range. However, when we calculate the susceptibility by 
adding a very small magnetic field term to the system 
(/i/Ji = 10^^), we find (Fig. |6|) that it diverges in the 
limit of zero temperature. To analyze this divergence, in 
the inset of Fig. we show the low-temperature behavior 
of x{T) on logarithmic axes. The data show a robust 
scaling relation of the form 



X' 



T- 



(8) 



where the exponent 7 = 1.0168(5). This scaling behavior 
indicates that the model is critical at zero temperature. 
Such zero-temperature critical behavior is also present in 
the antiferromagnetic q = 2 (Ising) pS], p9[ and q = 4 
Potts models on the triangular lattice pof, as well as the 
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FIG. 6: Susceptibility for the q — 4: Potts model on the gen- 
eralized decorated square lattice with Ji = J2 and D = 30. 
The inset shows the behavior near T = on logarithmic axes, 
which displays good linear scaling with a slope of 1.0168(5). 



(7 = 3 model on the square lattice M|. From our result 
we can infer that the critical point for the generalized 
decorated square lattice is not less than 4, i.e. Qc > 4. In 
comparison with the triangular lattice, where gc = 4 |30[| 
and the coordination z = 6, the lattice we consider here 
has a lower average coordination number, z = 16/3, but 
is also irregular in the sense of having inequivalent sites, 
an effect now understood [Q to increase qc. 

To conclude, we study the antiferromagnetic Potts 
model on the generalized decorated square lattice by 
calculating its field-dependent thermodynamic properties 
using the HOTRG method. For q = 3, we find a second- 
order phase transition, occurring at Tc ~ 0.528(1) for 
the isotropic point ( Ji = J2), to a state of partial order. 
The ground state has partial order on a single sublattice, 
with spontaneous breaking of 7r/2 rotational symmetry 
and the permutation symmetry of the Potts model. The 
entropy of this highly degenerate state scales with the lin- 
ear dimension of the system, meaning that the entropy 
per site vanishes at zero temperature. We obtain the 
phase diagram as a function of the ratio J2/J1, and our 
results indicate that a finite-temperature phase transi- 
tion is always present, but that Tc — > continuously as 
Ji -^ and discontinuously as J2 ^- 0. Thus the physics 
is determined not by the ratio J2/J1 but by the lattice 
connectivity. The same model with q = 4 is disordered 
at any finite temperature. However, for Ji = J2 we find 
that the model is critical at zero temperature, the sus- 



ceptibility diverging with the form x ~ T ^ , and thus we 
infer that the critical value of q in this lattice is qc > 4. 
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